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We examine a single-pulse preparation of the uniform su- 
perpositional wave function, which includes all basis states, in 
a spin quantum computer. The effective energy spectrum and 
the errors generated by this pulse are studied in detail. We 
show that, in spite of the finite width of the energy spectrum 
bands, amplitude and phase errors can be made reasonably 
small. 



1. Both the Shor and the Grover quantum algorithms 
begin with the preparation of a uniform superposition of 
the basis states. In Shor's algorithm, it is a superpo- 
sition in the x-register for the modular exponentiation: 
a^(mod7V). For Grover's algorithm, it is a superposi- 
tion of all possible entries of the unsorted data. In the 
language of computer science, the transformation of the 
ground state of the L-qubit register, 
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into the imiform superposition of all possible basic states, 

^unif = \nL~inL~2---nino), (nfc=0, 1), (2) 

is provided by the Hadamard transformation [Q . In phys- 
ical systems, this transformation can be implemented, 
for example, using two different methods: 1) by a se- 
lected 7r/2-pulse excitation of each qubit, and 2) by non- 
selective excitation of all qubits using a single 7r/2-pulse. 
The first method could be used, for example, for a chain 
of spins connected by the Ising interaction. Unfortu- 
nately, in this case non-resonant effects disturb the uni- 
form superposition. The reason is that a 7r/2-pulse acts 
not only on the chosen resonant spin but also on all other 
spins 1^. Besides, this method requires the application of 
L pulses. The second method requires only a single 7r/2- 
pulse. This method also cannot provide a perfect uniform 
superposition (2). The second method is used currently 
in developed statistical ensemble quantum computation 



In this paper, we analyze the second method. First we 
discuss the Hamiltonian of the system, then the effective 
energy spectrum in the rotating reference frame and, fi- 
nally, the error generated by a 7r/2-pulse. We present the 



results of numerical simulation of a chain which includes 
10 spins. We show that, in spite of the finite width of the 
energy spectrum bands, the amplitude and phase errors 
can be made acceptably small. 

2. Consider a chain of spin 1/2 nuclei described by the 
operators Ik- Assume that these spins have slightly dif- 
ferent Larmor frequencies, LOk and are connected by Ising 
interactions. In a liquid NMR quantum computation, 
one utilizes a statistical ensemble of such chains 
To prepare a uniform superposition (2), a 7r/2-pulse must 
be polarized along the (— ?/)-axis of the rotating reference 
frame (if initially the nuclear spins point in the positive 
z-direction) . The Hamiltonian of the system in the ro- 
tating frame is 0, 
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H = 5] [- K - + nil] -^JJ2 , (3) 
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where Ti — 1, is the amplitude of the pulse in the 
frequency units (the Rabi frequency), u> is the frequency 
of the pulse, and J is the constant of the Ising interaction. 
During the action of a 7r/2-pulse, the second term, 17/^, 
in the Hamiltonian (3) is the main one, the two other 
terms are supposed to provide small corrections to the 
spin dynamics. To choose the values of parameters, we 
assume that the following inequalities are satisfied for our 
spin quantum computer, 



(4) 



Assuming J/2tt ^ O.lkHz, (uik+i ~ ^^fc)/27r ~ IkHz, 
u!k/2TT ~ lOOMHz, we shall consider values of il/27r up 
to lOMHz. 

3. Next, we shall discuss the effective energy spectrum 
described by the Hamiltonian (3). To understand the 
behavior of the energy spectrum, we first analyze a sys- 
tem containing two spins only. When the inhomogeneity 
and the Ising interaction are absent, the energy spectrum 
consists of three lines: Eq = — O, Ei = 0, and E2 = n. 
The first level corresponds to the state, |00)_j,, where 
the index "-y" indicates that both spins point in the (- 
y)-direction ((— ?/)-representation). The twice degenerate 
level. El, corresponds to the states, |01)_j^ and |10)_y, 
and the energy level E2, refers to the state |ll)_j,. First, 
we consider the effect of inhomogeneity when the Ising 
interaction is absent. Assume, for example, that oj = luq 
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and u!i = luq + Auj. Then, the effective field for spin 
"1" in the frequency units is, Oe// = V^^^ + Aw^. The 
energy levels are given by the expressions, 

Eo = -{n + neff)/2, Ei^iQ- n,ff)/2, (6) 

E[ = {n.ff - n)/2, E2 = in + n,ff)/2. 

The main effect is the splitting of the central line into 
two lines with the energy separation, 

AE = E' -E^ {ALj)^/2n. 

Important fact is that this splitting decreases as l/ft as 
n increases. 

Let us consider the error generated by a non-resonant 
spin. If the initial state of this spin is \0)z, then with 
linear accuracy, Aw/fi, the wave function of the spin, 
'i'{t), can be written as, 

^{t) ^ [cos{nt/2)+i{ALj/n)sm{flt/2)]\0)z+sm{nt/2)\l)z 
After a 7r/2-pulse {ilt = tt/2), we have: 



and |10)_y, and E2 



-Eq. When the Rabi frequency. 
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Thus, in first order in Acj/il, the inhomogeneity gener- 
ates only a phase error which decreases as l/fi. 

Now, we consider the effect of the Ising interaction 
when the inhomogeneity is absent. At first sight, the 
Ising interaction produces an additional "effective field" 
in the z-direction, and its influence on the energy spec- 
trum must also decrease as fl increases. 



Q, increases, the influence of the Ising interaction on the 
energy levels, Eq and E2 decreases. But the splitting 
between the central energy levels does not change as 
increases: AEi = J. One might expect that this split- 
ting will generate an error which does not decrease as 
increases. Fortunately, it does not happen. If both 
spins point initially in the positive 2-direction, the wave 
function, ^'(i), can be written as. 
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+ze»^*/2(|01)_, + |10)_,)], 

where we neglected the terms ^ (J/fi)^. At the end of 
the 7r/2-pulse, we have. 
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FIG. 1. The energy spectrum as a function of fl. 

However, the "effective field" approach is not correct 
in this case. It is easy the check that the ground state is 
|00)_j^ with a small admixture of the state It has 

the energy, Eq = - v/f^^TTj/S^ • The central energy 
splits into two levels ± J/2, which correspond to symmet- 
ric and antisymmetric superpositions of the states |01)_jy 



Q. 

FIG. 2. The structure of a narrow strip of the central band 
shown in Fig|l[ 

For J = 0, we have a uniform superposition of the 
basis states (in both (— y)- and ^-representations) . For 
J ^ 0, to first order in J/f], the Ising interaction (similar 
to the inhomogeneity) generates only a phase error which 
decreases as l/fi. The non- vanishing bandwidth, AEi = 
J, does not cause a non-vanishing error. This happens 
because the phase error is proportional to the duration 
of a 7r/2-pulse which is proportional to l/fl. 

3. Next, we present the results of numerical simulations 
with L = 10 qubits. These simulations require operations 
on a digital computer in the Hilbert space with dimension 
D = 2^" = 1024. 
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FIG. 3. Widths AS6,(left) AE4, (right) of the 6th (central) 
and the 4th band vs Q. Open circles are for J = 0, crosses for 
J — 0.1. Dashed lines are the best fit to A/Q (circles) where 
A = 28.466 (left) and A = 23.673(right). 

If the Ising interaction and the inhomogeneity of the 
frequencies, Uk, are both absent, the energy spectrum 
of the Hamiltonian (3) consists of 11 equidistant levels 
separated by gaps with value Q. Both the interaction 
and the inhomogeneity cause splitting of all inner levels. 
This leads to formation of energy bands. 

Fig. 1^ shows the energy spectrum as a function of Q, 
for the following values of parameters, 

ujk+i - = 1, (fc = 0, 9), J = 0.1, ujk- LO ^ k - 4.5. 
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FIG. 4. Dependence on Q. of the maximal (open circles) 
and average errors (full circles ) for the amplitude modulus, 
|A„I, on n. The full line is the best fit, 0.2787/^^, while the 
dashed is 0.0953/Sl^. Other parameters are the same as in 
Fig.§. 

The last equation means that the frequency to of a 7r/2- 
pulse is equal to the average Larmor frequency, {uk) — 
(1/10) X]fe=o ^k- For each 17, there are 1024 energy levels. 
At the scale shown in Fig. |^, each band is represented 
by a point. 
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FIG. 5. Dependence of the maximal (open circles) and av- 
erage errors (full circles ) for the phase errors \cf>n\, on Q.. Full 
line is the best fit 13.0216/f2^hile dashed is 3.6606/f2. Other 
data are the same as in Fig.hl 

In fact, each band has a complicated structure. As an 
example. Fig. |^ shows the structure of a narrow strip of 
the central band as a function of VL. Fig. || (left) shows 
the dependence of the width of the central band, AiJg, 
on Vl (in the logarithmic scale) . When the interaction be- 
tween spins is absent (circles in Fig. ^ ), AE'g decreases 
approximately as \/VL. A finite interaction (crosses in 
Fig. 1^) changes this picture. After the width of the band 
caused by the inhomogeneity decreases to the value of ap- 
proximately 6 J (at K, 50), its value does not decrease. 
This dependence of AiJg on Q. is qualitatively similar to 
the results discussed above for two spins. Fig. ^ (right) 
shows similar dependence for the the width of the 4th 
band, A£'4. 
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FIG. 6. The distribution of the complex amplitudes. An, 
in the complex plane for Q, = 10'^ (left) and Vl — 10* (right). 
Other parameters are the same as in Figjl|. 

Now we consider the errors generated in the process 
of preparation of a uniform superposition of the basis 
states, 

In the rotating frame, the spin dynamics can be de- 
scribed as a superposition of stationary solutions with 
constant coefficients, which can be found from the initial 
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conditions. In the absence of both the interaction be- 
tween spins and the inhomogeneity, a 7r/2-pulse apphed 
to the ground state, |Ol_i...Oo)z, generates a uniform 
superposition of all 2^ basic states with the amplitudes 
l/%/2^. 
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FIG. 7. Level spacing distribution for a system of L = 12 
spins; Q = 100; u — loq, ujk = (^o + k, J = 20 (upper); J — 100 
(lower). The continuous line shows the Poisson distribution, 
while the dashed curve shows the GOE distribution. 

To describe the error in the complex amplitude. An = 
l^lnl exp(j$„), of the state |n)z (0 < n < 2^~^), we use 
two quantities: 1) 77 = |2~^/^ — |v4„|| which describes 
the error of the amplitude modulus, |A„|, and 2) the 
phase modulus, |<&n|, which describes the phase error (as 
= in the ideal case). Fig. ^ shows the dependences 
of the maximal error, r]max and the average error, r]ave, on 
the Rabi frequency, Q. One can see that both quantities 
decrease approximately as l/fi^. Fig. ^ shows similar 
dependences of the phase error on ft. Both, |$|maa: and 
l^laue decrease approximately as l/fl. This corresponds 
to the results derived above for two spins. Note that 
for parameters we chose (which correspond to J <C Aw), 
the main contribution to the error depends on the inho- 
mogeneity. The plots presented in Figs and || do not 
change significantly when J = 0. Fig. ^shows the dis- 
tribution of the complex amplitudes. An, in the complex 
plane for = 10^ and r2 = 10^. One can see that the 
distribution of An has a form of the arc whose length and 
width decrease as fl increases. 

4. In this section, we note that the quantum Hamilto- 
nian (3) belongs to the so-called quantum non-integrable 
systems. This means that, treated classically, these sys- 
tems exhibit chaotic behavior for some range of param- 
eters and initial conditions. (See, for example, |^ and 
references therein.) The quantum properties of classi- 
cally chaotic systems can be investigated by a detailed 
analysis of their eigenvalues and eigenfunctions |Q . This 
analysis will be the subject of future investigation. Here, 
we concentrate on the statistics of neighboring energy 
levels spacing, P{s). Indeed it has been conjectured that 
P{s) for a chaotic system depends only on some general 



symmetries g. In our case (for Hermitian matrices), we 
can assume the quantum chaos regime as described by a 
Wigner-Dyson distribution. Fig. ^ shows the transition 
of P(s), from the Poisson distribution to the Wigner- 
Dyson distribution of the Gaussian orthogonal ensemble 
(GOE) of Hermitian matrices. 

As one can see from Fig. |^, the transition to chaos ap- 
pears for relatively large values of J which are far beyond 
the range of quantum computation. Also, the appearance 
of the Poisson distribution, which should be a fingerprint 
of integrability, appears at relatively high values of J, 
out of quantum computation regime. Level statistics in 
the quantum computation regime will be the subject of 
future investigations. 

5. In conclusion, we investigated the errors generated 
by a single-pulse implementation of the Hadamard trans- 
formation for a chain of spins connected by the Ising in- 
teraction. In the rotating reference frame, the interaction 
between spins and the inhomogeneity of the Larmor fre- 
quencies split the energy levels into the bands. The char- 
acteristic width of the band caused by the inhomogeneity 
decreases as 1/ft as the Rabi frequency, fi, increases. The 
contribution to the band widths due to the Ising inter- 
action between spins remains constant. In spite of this 
fact, errors generated in the process of preparation of 
the uniform wave function decrease monotonically as 
increases. When the Rabi frequency, fl, increases, the 
errors of the amplitude modulus, \An\, decrease as l/fi^. 
The phase errors decrease only as l/fi. For reasonable 
values of SI, the errors of the amplitude's modulus, |A„|, 
become negligible. The phase error can be reduced to 
the order of 10~^rad. Thus, both the errors caused by 
the inhomogeneity and the interaction between spins can 
be made reasonably small, for a single-pulse generation 
of the superpositional wave function used in the main 
quantum algorithms. 
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